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We show that the ground state of a spin-1 Bose gas with 
an antiferromagnetic interaction is a fragmented condensate 
in uniform magnetic fields. The number fluctuations in each 
spin component change rapidly from being enormous (order 
N) to exceedingly small (order 1) as the magnetization of the 
system increases. A fragmented condensate can be turned 
into single condensate state by magnetic field gradients. The 
conditions for existence and method of detecting fragmented 
states are presented. 

Bose-Einstein condensation (BEC) was first intro- 
duced as a phenomenon in non- interacting Bose systems. 
The concept of BEC was generalized to interacting Bose 
systems in 1956 by Penrose and Onsager [Q. A system of 
N Boson is considered BE condensed if its single particle 
density matrix has one and only one macroscopic eigen- 
value (i.e. of order AT). The corresponding eigenfunction 
is identified as the quantum state macroscopically occu- 
pied. This characterization is in good agreement with 
the recent BEC experiments on magnetically trapped al- 
kali atoms, which are effectively scalar Bosons since their 
spins are frozen. Recently, optical trapping of Bose con- 
densate has become possible || . An optical trap confines 
all spin states, and the nature of the condensate depends 
on the magnetic interaction §{ B. In the case of 23 Na 
where the interaction is antiferromagnetic ||, the single 
condensate interpretation appears to agree with experi- 
ments ||. Since 23 Na is a spin-1 Boson, the macroscopi- 
cally occupied state is a three component spinor. 

However, in a recent paper, Law, Pu, and Bigelow ||] 
have pointed out that the Hamiltonian in ref. H {| in 
zero magnetic field has a singlet ground state, with prop- 
erties drastically different from those of spinor conden- 
sates. They, however, did not discuss how their results 
will reconcile with the MIT experiment |^] . The spin sin- 
glet turns out to be a "fragmented" condensate, mean- 
ing that its density matrix has more than one macro- 
scopic eigenvalue. The possibility of fragmented conden- 
sate was first discussed by Nozieres and Saint James (NS) 
[gj, who concluded that it cannot occur in homogenous 
scalar Bose gas with repulsive interactions. 

As we shall see, the singlet state is unstable. It 
quickly gives way to a much more generic fragmented 
state (called coherent-fragmented state) as magnetiza- 
tion increases. The coherent-fragmented state has iden- 
tical occupation number as the single condensate state in 
each spin component, except that there is no phase co- 
herence between different spin components. As a result, 
these two states cannot be distinguished by density mea- 
surements such as the MIT experiment M. To tell them 



apart, it is necessary to measure the phase coherence 
between different spin components, which can be done 
by performing Stern-Gerlach experiments along different 
axes. The origin of the fragmented state turns out to spin 
conservation. As a result, magnetic field gradients which 
destroy spin conservation can deform fragmented states 
toward single condensate states. The degree of deforma- 
tion depends on particle number and the strength of the 
field gradient. In the following, we first consider homoge- 
nous Bose gas, where "fragmentation" can be discussed 
most efficiently. Discussions on trapped gases will follow. 

I. Homogenous Spin-1 Bose Gas : Consider aspin- 
1 Bose gas with Hamiltonian ^^H = h + V, h = 
J>t/i(xW„ MxW = [-^V 2 + (7(x)]^- 7 B.S MI ,, 

^ = I / ^t^t^u^p [co^aii + c 2 S a/ 3 ■ Spv], where c , c 2 > 
0, M and 7 are the mass and gyromagnetic ratio of the 
Boson respectively, and B is a uniform magnetic field. 
The field operator Vv( x ) (/-* = =l=l 7 0) can be expanded 
as M (x)= £1~ 1//2 X^k^o eik Xfl M(k) where ^ i s the volume 
of the system. For simplicity, we shall denote a M (k = 0) 
simply as a^. Denoting the part of H containing a M alone 
as H and the rest as H ex , (H = H a + H ex ), we have 
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where S= a^S^ a u , N= a^a , and C= %[N 2 - AT]. 

To find the ground state \G), we first find the ground 
state of Hq (denoted as \F)) and then study conden- 
sate depletion effects due to H ex . From eq. ([!]), we see 
that \F) = \S total = S;Sl otal = S), where S is the in- 
teger closest to (jB£l/c2) — h, which is the minimum of 
(H ) F = (c 2 /2n)S(S + 1) - ~/BS - c 2 f + C. The equi- 
librium magnetization is s = S/N = jB/c2n + 0(A^ _1 ) 
and n = N/Q. In the following, we shall only consider 
the case 7B < c%n, where S/N ranges from to 1. 

In contrast, the optimum single condensate state in a 
magnetic field is 



\SC) 




vac) 
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where \ is a relative phase, N±\ = N(l ± y)/2, and 
y is the magnetization. (We shall from now on ab- 
sorb x m t° the operators for simplicity. The impor- 
tance of x will be discussed at the end.) The equilibrium 
magnetization y Q is obtained by minimizing the energy 
{H )sc = %N{N - l)y 2 - 1 BNy + C = G(y), and is 
y a = jBn/[c 2 (N - 1)] = jB/gn + 0(N~ 1 ). The energy 
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difference between \F) and |5C) is AE = (H ) SC -(H ) F 
= G{y )-G{s ) + {c 2 N/2n){2-s -s 2 ). Since s a < 1, the 
last term in AE is positive. Moreover, y — s Q ~ 0(A _1 ), 
wehaveG(2/ )-G( So ) = f^^-^ „ ^ wh i c h 
is smaller than the last term in AE by a factor of N and 
can therefore be ignored. This shows that for a homoge- 
nous Bose gas, the states \F) and |5C) are degenerate 
in the thermodynamic limit (A — > oo,f2 — > oo,N/Q — > 
finite), since their energies are of order A but their dif- 
ference (c2*n.) is of order 1. The relative stability between 
\F) over |5C) is therefore very delicate. To discuss this 
stability, it is necessary to understand the structure of 
\F), which turns out to be very remarkable. 

II. Super- and Coherent- Fragmentation : A sim- 
ple exercise shows that 0t = — 2a\o)_ 1 + aj 2 creates a sin- 
glet pair of spin-1 Bosons. The ground state \F) = |5; 5} 
is therefore given by 

\S;S) = -j=^= a\ s e^\ya,c), Q = (N - S)/2, (3) 
where f(Q; 5) is the normalization constant |?J 

fKni-wtf ^; 1 " . (4) 

Using eq.(^), it is easy to show that the single particle 
density matrix of \F) is diagonal, {p F ) a p — { a ^ a a )F = 
N a 5 a p, with 

A(5+l) + 5(5 + 2) (A-5)(5+l) 
1 25 + 3 ' 1 25 + 3 

Since p F has more than one macroscopic eigenvalue, \F) 
is a fragmented condensate for all 5 < N S. Again 
using eq.(|4|), the squared number fluctuation (AiVi) 2 = 
((a\a 1 — (aja^)) 2 ) for spin /i = 1 can be shown to be 

(ATVi) 2 = (53^3) (jj+k) + ((2I+3) 2 ) (jj+k) 

+ (^+5) ((fs+lf 2 ) • ( 6 ) 

Moreover, we have (AA_i) 2 = (AAi) 2 = (AA ) 2 /4 from 
the relations A_i= N% - 5, and A = A + 5 - 2N X . 

Eqs.(||) and (||) show that when 5 = 0, the system 
has N\ = No — N^i — A/3 Q with enormous fluc- 
tuations AN a ~ N. On the other hand, both Ao and 
AJV Q shrink rapidly as 5 increases. When 5 becomes 
macroscopic, No and AA Q become order 1, whereas N±± 
remains macroscopic, N±i — > (A ± 5)/2. The exceed- 
ingly small fluctuations AA Q means that state |5; 5) can 
be well approximated by 

|g;S)~pVi,^-i) = - 1 Ivac). (7) 

1^/ iVi!iV_i! 



To distinguish different fragmented states, we call those 
with (AA Q ) ~ A "super" -fragmented states, and those 
with (AN a ) ~ 1 "coherent" -fragmented state. In the 
thermodynamic limit, super- fragmented is a singularity 
which occurs only at S/N = 0. All other states with 
5/A ^ are coherent-fragmented. 

The origin of the small fluctuations in the coherent- 
fragmented state can be understood as follows. Not- 
ing that 5 + = \/2(aIao"'" a o a -i)' ec l-(0) can be written 
as H Q = H A +H B , 

Ha = ^ [5 2 + (2A - \)(N X + N_ ± )] - 7 S5 Z , (8) 

and Hb= ff( a o 2a i a -i +h.c). Hb contains all terms 
which are responsible for transformation among spin 
species. In Ha we have dropped terms depending only 
on A. Clearly, Ha is minimized when Ao = 0. If 
I Ai, A_i, No) denotes the state with A M Bosons with spin 
/i, the ground state of H A is |A 1 ,A_ 1) 0)= \N 1 ,N_ 1 ), 
where 2V±i = (A ± 5)/2, and 5 = jBfl/c 2 . The ground 
state energy is E = — c 2 5 2 /(2fi). The effect of Hb is 
to mix in N a ^ states \q) = \N 1 - q,N_ 1 - q,2q). To 
leading order in N, we have H Q — E 1 = H, 

H= C ^~ £ [2q\q){q\ + \ q (\q + l){q\ + h.c.)} (9) 



where X q ^(Ai - q)(N_ 1 - q)(2q + 2)(2q + Tj/N For 
1 « q « A 1; A_!, A,= x(q+l), x =y/l - (5/A) 2 . To 
estimate the energy of eq.(pl), we consider the following 
variational state |\I>) = y/a l) q e~ aq / 2 with energy 

E =^^-{^—xe~ a / 2 [—+§])■ The minimum condition is 
e a/2 = x ^ _|_ Iq,^!^^ rp Q j ea( j m g 0( j er j n 5*/tv, we have 

a = v2(5/A). The energy correction (£T)* is therefore 
of order ^q-, which is lower than E a by a factor of A. 
Thus, to the leading order in A, H D can be replaced by 
Ha [eq.d|)] with N = 0, with ground state eq.(^). 

For later discussions, it is useful to compare 
the coherent-fragmented state eq.(fTj) with the sin- 
gle condensate state \SC) eq.@. Writing |5C)= 

E^-iv, ( iSwfe^yi ^j ^ where ^ = ^ + 

£, N_i — £), and using the Stirling formula, it is straight- 
forward to show that the single condensate is a Gaussian 
sum of coherent fragmented states 

N-i 

\SC) = (na 2 )- 1 /* e- e2 / 2 " 2 e^\e) (10) 

e=-Ni 

where a 2 = 2A^iA_i/A, £ = (A^ 1 - Af^/4 [|. Within 
the space of /x = ±1, the density matrix of \SC) is 
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The off-diagonal elements y/NiN-i are absent in p F . 

III. Persistence of Fragmentation : So far, we 
have ignored condensate depletion (i.e. due to H ex ), 
and various realistic atomic physics effects, (see below). 
In the following, we shall focus on coherent-fragmented 
states as they are most common. Following the method 
of Huang and Yang |10| to extract the dominate con- 
densate depletion effect from H ex , we find that the effec- 
tive Hamiltonian for condensate depletion for both single 
condensate and fragmented state have identical Bogoli- 
ubov form, which is obtained by replacing the operators 
Gy-i(k = 0) by the c-number y/N±i- This means that 
the thermodynamic degeneracy of \F) and \SC) cannot 
be lifted by condensate depletion, as it gives rise to same 
energy change in both cases. (This result also applies to 
the trapped Bosons, in which case the states (k, — k) are 
replaced by opposite angular momentum states.) As for 
atomic physics effects such as mixing of different hyper- 
fine states and quadratic Zeeman effects, they all respect 
spin rotational symmetry along z [[TT|| . The new ground 
state \G) with all these effects included therefore remains 
an eigenstate of S z , which implies (G\a\a_ 1 \G) = as 

a\a_ 1 changes si°\ Thus, these effects cannot reassem- 
ble a fragmented state into a single condensate state, for 
(a\a_ 1 )sc — -JN1N-1 is of order ~ N instead of 0. 

IV. Reassembling Fragmented States into Sin- 
gle Condensates : Since spin conservation protects 
fragmentation, we consider perturbations that break spin 
symmetry. The natural candidate is a magnetic field 
gradient G' . To be concrete, we take B(x) = B Q (z + 
G'[x3c — zz]). The condition of small field gradient is that 
G'Sl -1 / 3 << 1. It is useful to choose local field direction 
B(r) as the spin quantization axis. This is done by per- 
forming a unitary transformation U = YiiLi e~^ Xi )' s % 
where = z x B = G'xy + 0(G' 2 ). The interaction 
V is invariant under U because it is a spin conserving 
contact interaction. However, WpU = p + hG'S y x + 
0(G' 3 ). This causes H a -> (U^H U) = H Q + H u 



where H\= eX)»=i ) 2 anc ^ e = h 2M ■ When oper- 
ated on the coherent fragmented states eq.(^), H\ is re- 
duced to Hi = -(e/2)(a[a_ 1 + al^) + eN/2. The ef- 
fect of Hi on the state |JVi,JV_i) is to generate the set 
{\£) = \N t + £,N_ X - £)}. Within this set, H' has the 
tight-binding form 



H'o 



E 

=0,±1, 



c ^\m-\{\i- 



!)(£{ + h.c. 



(12) 



where t e = e y /(N 1 + £ + l)(A_i - I). The eigenstates of 
eq.(|l^), I*) = Yle^Afyi satisfies the Schrodinger equa- 
tion E^ e = (2c 2 /n)e 2 V e - (ti9 t+1 + t<_i*i_i)/2. Al- 
though this equation can be solved numerically, it is more 
illuminating to consider the following analytic approxi- 
mation, which turns out to be very accurate. As we shall 



verify later, the number of I terms in the ground state 
|\E r ) is much less than the typical value of N\ and N—\, 
so that we can replace (te + te-i)/2 ~eWNiN-x, and 

(te - t t -i)/2 e^NTi^ where £ = (AC^ - N^)/A 

as defined before. The Schrodinger equation in the con- 



tinuum limit is then (j^+r] 4 )^e= 



14 d^ e , ^4 
2 'I 



dl 2 



where 



V 4 = e^N^Cln/^cz). 
The normalized ground state is 



i*> = (-7 2 r 1/4 £' 



-t*/2r, 2 M 



\£), 



(13) 



(14) 



with a density matrix (p^r) a /3 = (^\ a a a g\^) ( a iP = 



(Maj3 



Ni 



v /AW~i"e~ 1/4 ' ? 
The eigenvalues of p^ are 
1 



A+ = - 



N ± yj N^M 1 + S 2 (l - ^W) 



(16) 



For zero field gradient, r\ — > 0, \± — ► ^(N±S). For large 
field gradients, rj » 1, eq.Q) reduces to eq.(||), and 
A+ — > N, and A_ — * 0. The systems turns into a single 
condensate. Note that even for r\ ~ 5, the system is es- 
sentially a single condensate state. Using the expression 
of C2 in ref. ||, C2 = AnTi 2 Aa sc /M , Aa sc — (02 — ao)/3, 
we have 77 =r] (l - (^)) 1/8 , and ?? = [eftA/(4c 2 )] 1/4 

Since 77 ~ A 1 / 3 , arbitrarily small 



(g'st 



N 



1/4 



327r(Ao sc n 1 / 3 

field gradients will change a fragmented state into a sin- 
gle condensate for homogenous Bose gas in the thermo- 
dynamic limit. 

However, trapped Bose gases with N ~ 10 are 
in mesoscopic rather than thermodynamic limit. In 
this limit, super- and coherent-fragmented states are no 
longer singularities in the phase space (jr,e). Instead, 
they occupy a finite region in the phase space which 
crosses over to single condensate states in a continuous 
manner. This implies the possibility of observing frag- 
mented condensates in trapped gases. 

V. Trapped Spin-1 Bose Gas : When B ^ 0, differ- 
ent spin components have different spatial extents. For 
cylindrical traps, we write Vv( x ) = /^( x ) a ^i + <?V( X )> 
where / 's are normalized and cylindrically symmetric 
wavefunction of the condensate with spin //, and is the 
non-condensate part of ip^ satisfing J f fl cf) = 0. The exact 
forms of fn are determined by energy minimization. As in 
the homogenous case, we write H = H Q + H ex where H Q 
contains only a^. H D can again be written as Ha+Hb, 
where Ha contains only S z , N t + N_ 1 , N Q ; and Hb has 
the same form as before but with a different coefficient. 
Ha is again minimized at No = 0. The effect of Hb is 
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again to give rise to eq.([|) with a weaker hopping term. 
||12fl Our previous results therefore apply, i.e. to lead- 
ing order in N, Hb can be ignored and the ground state 
with macroscopic magnetization is accurately given by 
|F)= (iV^ALi), N ±i = (N±S)/2. The functions f x and 
/_i are determined by minimizing E — (H )p. In the 
presence of a field gradient, a repeat of our previous cal- 
culation shows that the Hamiltonian H' a is again given by 

eq.d), with tt" 1 - i / [(A 2 + f 2 ,) 2 + ^(f 2 - f\f 
and e -> ef 

To estimate the field gradients G' and magnetization 
S/N needed for observing fragmented states, we recall 
from eq.(||) that t? 4 = 3%^" (1 - (#) 2 ) 1/2 , where we 
have used C2 = Airh 2 Aa sc /M . Next, we note that the 
maximum and minimum value of A+ in eq. (|l6|) is N and 
(N + S)/2. Taking the mean of these values A* = (N + 
(N + S)/2)/2 as a specification of the crossover region 
between single condensate and fragmented state, we find 
that fragmented states with A+ < A* emerges when G < 

^ \ / M£+fiy For a 23 Na gas 

(Aa sc ~ 10~ 8 cm) with N = 10 6 , S/N = 0.2, we have 
A* /N — 0.8. Since f£ is roughly the volume of the system, 
if we take for £1 / ~ 10 _12 cm 3 , we have A < A* 

when G" < 0.46cm- 1 ) ||. 

Finally, we note that coherent-fragmented states cross 
over to super-fragmented states as S decreases. Set- 
ting AN a ~ y/~N as the crossover estimate from super- 
to coherent- fragmentation. Eq.(||) implies that (for 
N, S >> 1) super-fragmentation emerges when S/N < 
1/V8N. For N = 1800, it means S/N < 0.016. Thus, 
super-fragmentation can only be achieved when S/N is 
very close to zero. On the other hand, coherent frag- 
mented states are more generic and are easier to realize. 

VI. Distinguishing fragmented states from sin- 
gle condensate states: Because of the difference 
in density matrices, the states \SC) and \F)= \S;S) 
have different rotational properties. Under a spin ro- 
tation (with Euler angles (a, (3, 7)), the density matrix 
Piiv = { a t a fj) becomes DpD\ where D{a,j3,^f) is the ro- 
tational matrix. The number of spin-/x particles changes 
from — p^f! to N' — {DpD^)^^. For single conden- 
sate state of the form eq. (Q) before rotation, we have 

N ±1 = c 4 N ±1 + s 4 N Tl + 2c 2 s 2 v / iViiV-iCOs2x (17) 



Nn 



2s 2 c 2 (N 1 + N_{) - 4c 2 s 2 v /7Vi7V_iCos2x (18) 



where \— 7 — x/2, s — sin/3, c = cos/3, and (3 is the angle 
between the original and the new quantization axis. 

For fragmented states with occupations numbers 
N\ , No , iV_ 1 before rotation, we have 



eqs.(|17]) to ( |20| ) show that the difference in N' between 
single condensate state and fragmented state is the phase 
coherence term cx sin 2 /3-\//Vi/V_iCOs2;\\ This difference 
can be detected by first preparing a sample with finite 
magnetization along a specified axis and then performing 
a sequence of Stern-Gerlach experiments along a different 
axis, as well as different axes. 
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2 1 / 2 ) 5 (d x + id y ) s (V 2 ) Q , and 
vac) 

', we have f(Q, S)= {Q, S\Q, S} = [T* c3 (x)T«(x') 



Ay 
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These statements can be seen most clearly in 
the local density / Thomas-Fermi approximation. 
The ratio between \ q and the diagonal "potential" 
term is given by, for 1 < < q < < A_ 1 , Ai ; 
v^^vT / / 2 /i/-i / / / 2 [iVi/i 2 + A-i/^x] < 1. 
The effect of quadratic Zeeman effect on the crossover 
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N = 2s 2 c 2 (N 1 + JV_i) + A (c 2 - s 2 ) 2 
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In the coherent-fragmented regime where N±i >> Nq, 
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